function sensitivity_analysis(base_params)
    % 参数敏感性分析
    parameters = {'A', 'sigma2', 'L', 'T'};
    variations = 0.8:0.1:1.2;  % ±20%变化
    
    figure('Position', [100, 100, 1000, 800]);
    for i = 1:length(parameters)
        subplot(2,2,i);
        legend_entries = cell(length(variations), 1);
        
        for j = 1:length(variations)
            v = variations(j);
            params = base_params;
            params.(parameters{i}) = params.(parameters{i}) * v;
            
            % 计算散射结果
            [theta2_range, ~, yy] = calculate_scattering(params);
            
            % 绘制结果
            plot(rad2deg(theta2_range), yy, 'LineWidth', 1);
            hold on;
            
            % 准备图例
            legend_entries{j} = sprintf('%.1fx', v);
        end
        
        % 设置子图属性
        title([parameters{i} ' 变化的影响']);
        xlabel('散射角度 (°)');
        ylabel('log_{10}|ρρ^*|');
        grid on;
        legend(legend_entries);
    end
end 